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Abstract. The existence of positive equilibrium solutions to age-dependent population equations with non- 
linear diffusion is studied in an abstract setting. By introducing a bifurcation parameter measuring the intensity 
of the fertility it is shown that a branch of (positive) equilibria bifurcates from the trivial equilibrium. In some 
cases the direction of bifurcation is analyzed. 



1. Introduction 



The present paper is dedicated to the study of nontrivial equiUbrium (i.e. nonzero time-independent) so- 
lutions to abstract age-structured population models with nonlinear diffusion, that is, to equations of the 
form 



subject to some initial condition at i = 0. Here, u = u{t, a) is a function taking on values in some Banach 
space Eq and represents in applications the density at time t of a population of individuals structured by 
age a £ J [0, a„i), where a,n £ (0, oo] is the maximal age. The real-valued functions /i = ii{u, a) 
and (3 = P{u, a) are respectively the death and birth modulus. The operator A{u, a) depending in a certain 
way on the density u specified later governs the spatial movement of individuals. It is assumed to be a 
(unbounded) linear operator A{u, a) : Ei C Eo ^ Eq satisfying additional technical assumptions given 
later. 

Age-structured models have a long history and questions regarding well-posedness and behavior for 
large times were investigated (see ||3T1 and the references therein) though most research was devoted to 
models neglecting spatial structure from the outset or considering merely linear diffusion, see e.g. lfT2l[T4l 
[161 [TSl l22l |26l and the references therein. Less seems to be known for the case of age-structured models 
with nonlinear diffusion (however, see e.g. P1 [T5ll28ll29l[T7l ). 

The understanding of the large time behavior of age-structured populations whose evolution is governed 
by equations ( 1 1.1b . (I1.2l i requires in particular precise information about the existence of equilibrium solu- 
tions. Since obviously m = is such an equilibrium solution the aim is to establish existence of nontrivial 
equilibria. Moreover, since u represents a density the main task is to single out the positive equilibrium 
solutions in the (ordered) space Eq. 

Equilibria of (ll.lt . jl.lj are solutions to 



dtu + daU + A{u, a)u + a)u = , 



t > , a e (0, am) , 



(1.1) 




t>0, 



(1.2) 



daU + A{u, a) u + fi{u, a) u = 



a e J 



(1.3) 




(1.4) 
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Supposing that the map (with u fixed) 

a i-^ A{u, a) :— A{u, a) + fi{u^ a) 

generates a parabolic evolution operator lit, (a, (y),Q < a < a < a„i, on Eq, a solution to (11.31 1. ( 11.4b must 
satisfy 

u{a) = n„(a, 0)?i(0) , a € J , u{0) = Q{u)u{0) , (1.5) 

that is, u(0) is (if nonzero) an eigenvector corresponding to the eigenvalue 1 of the linear operator Q{u) on 
Eq (for u fixed) defined by 

Q{u):^[ I3{u,a)nu{a,0)da . (1.6) 



^0 

Roughly speaking, Q{u) contains information about the spatial distribution of the expected number of 
newborns that an individual produces over its lifetime when the population's distribution is u. In the present 
paper we suggest a bifurcation problem by introducing a bifurcation parameter n which determines the 
intensity of the fertility without changing its structure. More precisely, we are interested in nontrivial 
solutions {n, u) (that is, m ^ 0) to 

daU + A(u, a)u + ^{u, a)u ~ , a ^ J , (1-7) 
u(0) — n b{u,a)u{a)da , (1-8) 







where we put 

nb{u,a) := f3{u,a) , (1-9) 
with b being a normalized function such that the spectral radius of the bounded hnear operator 



Qa := / 6(0,a)no(a,0)da 
Jo 

equals 1, that is, 

r(Qo) = l. (1.10) 

Note that under this normalization we have r{Q{u)) ~ nr{Qo) — n; the bifurcation parameter n is thus 
the spectral radius of the "inherent spatial net production rate at low densities" (technically when u = 0). 
If r{Qo) is an eigenvalue of Qo, then dl.lOt may be interpreted as that there exists a distribution for which 
the population is at balance meaning that the birth processes yield exact replacement (provided that death 
and birth modulus and spatial displacement are described by /i(0, •), /3(0, •), and A{0, •)). 

In this paper we provide a set of n values for which (11.7b . dl.Sb have nontrivial and positive solutions, 
respectively, around the critical value n = 1 and u = 0, analogously to the "spatially homogeneous" case 
(i.e. when A = 0), see |6 |. More precisely, it is shown that a branch of nontrivial solutions bifurcates from 
(i.e. intersects with) the branch of trivial solutions (n, u) = {n, 0), n G R, at the critical value of n. In 
principle, the direction at which bifurcation occurs will be related to (the values at li = of the derivatives 
of) /i, (3, and A computing a parametrization of the branch of nontrivial solutions. In some cases, the 
direction can be computed explicitly. In particular, examples will be given where supercritical bifurcation 
occurs. 

In order to derive the results we consider problem ( 11. 7b . dl.Sb in a more general framework so that the 
results actually apply to a broader range of similar problems. In Section|2]we investigate the general abstract 
framework and prove the bifurcation result using findings based on the implicit function theorem obtained 
in Q. We also derive a more precise characterization of the nontrivial branch of solutions and show that the 
equilibria are positive. The subsequent Section [3] then gives applications for these results. We shall point 
out that analogue results for populations structured by age only (that is, when A — 0) were derived |l6][7][8]. 
Furthermore, additional results regarding equilibrium solutions for age-structured equations are to be found 
in e.g. ||T0llin[I3[Il[l9l|20l|2T]|30l|3l] and the references therein. 
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2. Abstract Formulation 

Given Banach spaces E and F we write C{E, F) for the space of bounded linear operators from E to 
F equipped with the usual operator norm, and we put C{E) := C{E, E). We write r{A) for the spectral 
radius of A e ^{E). The subspace of C{E) consisting of compact operators is IC{E). If E is an ordered 
Banach space we write respectively C+{E) and 1C^{E) for the corresponding positive operators. We let 

Cis{E, F) denote the subspace of C{E, F) of topological linear isomorphisms. If £' > F, that is, if E 
is densely embedded in F, then T-L{E, F) is the set of all negative generators of analytic semigroups on F 
with domain E. BIV{E; <f>) stands for the set of operators with bounded imaginary powers and power angle 
4> e [0,7r/2), that is, those linear operators A in i? for which there is A/ > 1 such that ||^'*||£(£;) < Me"^!*!, 
t e M. For details we refer to Q. 

Throughout this paper we suppose that Eq is a real Banach space and Ei Eq, that is, Ei is a 
densely and compactly embedded subspace of Eq. We fix p e (l,oo), put c, c,[p) := 1 — 1/p and 
set E^ := {Eo,Ei)^^p with (•, •)^.p being the real interpolation functor Similarly we choose for each 
a G (0, 1) \ {1 — 1/p} an arbitrary admissible interpolation functor (•, •)„ and put Ea {Eq, Ei)a so 

d , 

that El Ea (see [2J). If Eq is ordered by a closed convex cone Eq, then the interpolation spaces are 
equipped with the order naturally induced by E^ . Given G (0, oo] we set J := [0, a^) which thus may 
be bounded or unbounded. Moreover, we put 

Eo LpiJ, Eq) , El := Lp{J, E^) n W^J, Eq) 

and recall that 

El BC/C(J,F,) (2.1) 
according to, e.g. lH III.Thm.4.10.2], where BUC stands for bounded and uniformly continuous. In par- 
ticular, the trace := u(0) is well-defined for u € Ei. 

We then study problems of the form 

daU + A(?i, a) u = , a € J , (2.2) 

u(0) = ni{u) (2.3) 

where A{u, a) e C{Ei, Eq) and £{u) E Eq for a £ J and u G Ei with l{{)) = 0. We will impose more 
restrictions later Introducing Ao(a) :— A(0, a) and assuming a decomposition 

e{u) = iQiu) + e,{u) 

with a linear part Iq, we first focus our attention on the linearization around of the above problem. 
2.1. Preliminaries. In the following we assume that 

£q e C{Ei,Ei)) for some d E (<r, l) (2.4) 

and that 

Ao E LooiJ, jO.{Ei, Eq)) generates a paraboUc evolution operator 

Ho (a, fj), < cr < a < on Eq with regularity subspace Ei and (2.5) 
possesses maximal Lp-regularity, that is, {da + Aq, 7) G £is(Ei, Eq x E^) . 
For details about evolution operators and operators possessing maximal regularity we refer the reader, e.g., 
to It seems to be worthwhile to point out that, owing to ( 12. 5t and |2 III. Prop. 1.3.1], the problem 

daU + Ao(a)u = /(a) , a E J , u{0) — 
admits for each datum (/, E Eq x E^ a unique solution u G Ei given by 

u{a) = UQ{a,0)u° + [ UQ{a,a)f{a) da , aEJ, (2.6) 
Jo 
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satisfying for some cq > 

\\u\\E,<Co{\\f\\Eo + \\u°\\E,) . (2.7) 

In particular, 

Uo{-,0)eC{E„Ei) and e /:(IEo, Ei) , (2.8) 

where 

{Kof){a) := / Ho (a, a)/(a) da, a€J, / £ Eq , 

Jo 

while, due to (|2.4| i. 

(/^4(ifo/)) e/:(Eo,£^.) . (2.9) 
Moreover, we obtain from ( l24l i. (IZST i, and the fact that =^ (e.g. see m I.Thm.2. 11.1]) 

Qo e C{E,,E^)r]IC{E,) for Qow := 4(no(-, 0)w) , w G . (2.10) 

The next result will be fundamental for what follows. 

Lemma 2.1. Suppose ( 12.41 ) and i2.5\ . Then the operator 

Lu := {-fu - {da + Ao)u) 

satisfies L S £(IEi , iS^ x Eq) and has a closed kernel ker(L) and a closed range rg(L) of finite dimension 
and codimension, respectively, both of which admit bounded projections Pk and Pr- In fact, 

dim(ker(L)) = codim(rg(L)) = dim(ker(l — Qo)) < oo 

and 

ker(L) = span{no(-, 0)w ; w e kcr(l — Qo)} , 
rg(L)^ := (1 - Pr){E, x Eq) = rg(l - Qo)^ x {0} , 
where E, = rg(l - Qo) ® rg(l - Qo)^- 

Proof. First observe that ( 12.6b implies that, for {hi, 112) £ E^ x Eo, the equation Lu — (/ii,ft.2) with 
M G El is equivalent to 

u = no{;0)u{0) + Koh2, {1 - Qo)u{0) = hi + eo{Koh2) . (2.11) 

If 1 is not an eigenvalue of Qo G K.{E^), then ( 12.1 It easily entails that kcr(L) is trivial. Moreover, in this 
case we have hi + ^^0(^0^2) G E^ for any {hi, /i2) G -E^ x Eo by ( 12.91 ) and there is a unique w E E^ for 
which (1 — Qo)w = hi+ £o(-f^o^2)- Thus u := IIq{-, 0)w + Koh2 belongs to Ei due to ( 12. 8t and satisfies 
Lu ~ {hi, h2), whence rg{L) — E^ xEq from which the claim follows in this case. 
Otherwise, if 1 is an eigenvalue of Qo G IC{E^), then (12.1 lb ensures 

ker(L) = span{no(-, 0)w ; w G ker(l - Qo)} C Ei 

which is clearly closed since L G >C(Ei, E^ x Eo) by ( 12.41 ), ( 12. 5b . In particular, the dimension of ker(L) 
equals the dimension of ker(l — Qo), the latter clearly being finite since the eigenvalue 1 has finite mul- 
tiplicity. Therefore, ker(L) is complemented in Ei and admits a bounded projection Pk G C{Ei , ker(L)) 
(e.g., see {23^ Lem.4.21]). Next, given {hi,h2) G rg(L) C E^ x Eo and u G Ei with Lu = {hi,h2), 
we have hi + (-Q{Kah2) G rg(l — Qo) as observed in ( 12.111 ). Conversely, if {hi, ^,2) G E^ x Eo and 
(1 — Qo)w = /ii + lo{KQh2) for some w G E^^, then Lu — {hi, h2) for u :— I\-o{-, 0)w + KQh2- Thus 

rg(L) = {{hi,h2) G E, X Eo ; hi + 4(^^0/12) G rg(l - Qo)} • (2.12) 

Since Qo is compact, M := rg(l — Qo) is a closed subspace of E^, hence rg{L) is closed in E^ x Eo due 
to ( 12. 9b . Furthermore, codim(M) = dim(ker(l — Qo)) < 00, and hence M is complemented in E^, that 
is, E^ = M © A/-'-. Let Pm e C{E^) denote the projection onto M along A/-'- and set 

Pr{hi,h2) := {PMhi - {l-PM)eo{Koh2),h2) . (2.13) 
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Then clearly ^ Pr e C{E, x Eq) by ^M, Pr{E, x Eq) = rg(i) by jTHI l. and 

(1 - P.)(/ii, /12) = ((1 - PM)ihi + £o{Koh2)),0) e X {0} . 

Thus we conclude that i?^ x Eq = rg{L) ® Tg{L)-^ with rg(L)^ ~ x {0}, and so 
codim(rg(i)) = dim(M-L) = dim(kcr(l - Qo)) = dim(kcr(L)) . 

This proves the assertion. □ 

The verification of i2.5i is not a simple task in general. We thus recall conditions that allow us in 
Section |3]to consider cases for which i2.5i is readily verified. 

Lemma 2.2. Suppose that 

Aq G BUC{J, C{Ei, Eq)) generates a parabolic evolution operator on Eq . (2.14) 

Further suppose, for each a G J, that belongs to the resolvent set o/A(a), that 

Ao(a) possesses maximal Lp-regularity , (2.15) 

and that 

lira Ao(a) exists in C{Ei, Eq) if am = 00 . (2.16) 

a — >oo 

Then ( 12.5b is satisfied. 

Proof. This is a consequence of 1241 Thm. 1 .4] . □ 

Remarks 2.3. (a) If Aq £ CP{J,Ti.{Ei, Eq)) for some p > 0, then it generates a parabolic evolution 
operator on Eq due to ||2l II.Cor.4.4.2]. 

(b) In case that Eq is a UMD-space (see |2| for a definition and properties), condition (12. 15b holds if for 
each a £ J there is some angle 6{a) G [0, Tr/2)forwhich Ao(a) G BTP{Eq; 0(a)), see [2, III. Thm.4. 10.7]. 

2.2. Nonlinear Theory. We now focus on problem ( 12. 2b . ( 12. 3b . Let m G N \ {0} and let E denote an open 
ball in Ei centered at of some positive radius Rq > 0. Suppose that 

AeC"^{^,L^(J,C{Ei,EQ)) and Aq := A(0) satisfies (ES ■ (2.17) 

We set A, (w) := A(u) — Aq and sometimes write A(u, a) :— A(u)(a) for u G S, a G J and accordingly 
A,(m, a) := A*(u)(a). We also assume that £ admits a decomposition 

e(u) = £q(u) + £4u) , (2.18) 

where the hnearpart £q satisfies ( |2.4b and £^ is such that £^(eu) = e£^,(e,u), w G S, |e| < 1, for some 
function 

G C7™((-1,1) x !],£;,) with Z»(0,-)-0, Dj.,(er)^0- (2.19) 

We put 

T(\,u) X{£q(u),0) + ((A + l)4(w),-A,(ii)ii) , (A, u) G M x S , 

and note that with n = A + 1 problem (12. 2b . (12. 3b can be be re-written as Lu = T(A, u) with L being given 
in LemmaO We then introduce T G C"'(M x (-1, 1) x S, E, x Eq) as 

r(A, £, u) := X{£q{u),0) + ((A + l)£,(e, u), —A^,{eu)uj , (A, e, u) G R X (—1, 1) x S , 

and observe that T(X,eu) — eT(X,e,u). Nontrivial solutions to (12. 2b . (12. 3b are then provided by the 
following 
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Theorem 2.4. Suppose ( 12. 4t , ( |2.17t , ( 12.181 ), ant/ ( |2.19t . Moreover, suppose that r{Qo) = 1 is an eigen- 
value of Qq G lC{Et;) with geometric multiplicity 1, where Qq is defined in ( 12.1 Oi l, ant/ let B be a 
corresponding eigenvector Then there exists Eq > such that the problem 

daU + A(u, a)u — , a ^ J , 

u(0) = ne{u) 

has a branch of nontrivial solutions | {n{£), u(e)) e R'^ x Ei ; < |e| < eo} of the form 

n(e) = l + A(e), u(e) = e(no(-, 0)5 + , < |£| < eq , 

where A : (— eoi^o) ^ ^ cind z : (— eo,£o) ^ ker(L)^ are m-times continuously differentiable with 
A(0) = Gflnc/z(0) = 0. 

Proof. We re-write ( |2.2| i. ( 12.31 ) as Lu = T(A, u) and validate the requirements for Theorem 1 in ||5]- First 
recall that Lemma IZTI warrants that L e £(Ei, x Eq) has a closed kernel ker(_L) = span{no(-, 0)5} 
and a closed range rg(L) both admitting bounded projections Pk and P^, respectively, and that the codi- 
mension of rg(L) equals 1. Thus HI and H2 in [5] hold. To validate H3 therein we just have to observe 
that for y e ker(L) n S 

T(0,0,2/) = (4(0,j/),-A,(G)y) = (0,0) 

and 

D3T(0,0,z/)= (0,-A40)y) = (0,0). 
It remains to verify H4 in jS). For, let 1 — be the projection of Ec; x Eq onto the one-dimensional space 
rg(L)-L = X {0} with M = rg(l - Qo) and let c(A, e, z) be the component of T(A, e, no(-, 0)5 + z) 
with respect to the basis {(P, 0)} of rg(L)^ for given A G M, |e| < 1, and ||z||ei < ^o/2. Here we may 
assume without loss of generality that ||no(-, 0)i3||Ei < ^o/2- Hence it follows from Q^B = B E M'^, 
(irT9] l. and dTlJl l that 

(i-p,)r(A,o,no(-,o)p) = (i-p,)(AB,o) = A(B,o) , 

that is, c\{X, 0, 0) = 1. Now ^ Thm.l] implies the assertion. □ 
Remark 2.5. Clearly, the result applies to non-homogeneous problems 

daU + A{u, a)u = g{u, a) , a J , 
u{0) = n£{u) 

as well provided that there is g such that g{eu, •) — sg{e, u, •) with 

[(e,u) H^5(e,u,-)] G (7™((-l,l) X S,Eo) and 5(0,-,-) = 0. 

Next we compute the e-expansion of the branch (n(e), u(e)). Under the assumptions of Theorem |2.4| we 
write El = ker(L) © ker(i)^ and let Pk G >C(Ei) denote the projection onto ker(L) = span{no(-, 0)P} 
such that PfcU = /t(-u)no(-, 0)i? with l^{u) G M for u G Ei. Again we set M = rg(l - Qq) and 
E^ ^ M ® with corresponding projection Pm (see the proof of Lemma l2.1b . 

Proposition 2.6. In addition to the assumptions of Theorem \2.4\ with m > 2 suppose that D^i* (0) = 0. 
Then the branch of nontrivial solutions {n{e), u{e)), \e\ < eo,from Theorem \2.4\ can be written in the form 

n{£) ^l + Ce + n^e) , ^(e) = £no(-, 0)P + (no(-, 0)^ - i^o/^) + ew,(£) 

for \e\ < Eq, where n^: : (— £o,£o) ^ and u^, : (— eo,eo) ]ieT{L)^ are such that \n^,{e)\ — o{e^) and 
||w*(e)||Ei — o(£^) as \e\ — > 0. The function /i G Eg is definedby 

h{a) := (i?„A,(0)(no(-,0)P)(a))no(a,0)P, aG J, 

C G M the unique coefficient of 

(1 - PM){io{Koh) - g) ^(Be 
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with 



9 lDle4Q)[Uo{;0)B,Uoi;0)B] e E, , 



and £ Ei; is the unique solution to 

(1 - Qo)e = CB + g- eoiKoh) e M , ^(no(-, 0)0 = tiKoh) . 

Proof. We plug the twice continuously differentiable functions A = A(e) and u = u{e) provided by Theo- 
rem l2.4l into the equation Lu = T(A, u) which we then differentiate twice with respect to e. Evaluating the 
result at e = and using D„£*(0) = together with £o(^o{', 0)B) ~ B, we obtain 

Lz'(O) = (A'(O)B + g, -h) (2.20) 

with dashes denoting derivatives with respect to e and g,h as given in the statement. Hence, from (12.1 II) . 

y:^ X'iO)B + g-eoiKoh)eM 

and thus, since PmB = 0, 

{1 - Pm){ - g + io{Koh)) = \'{0)B 

from which the formula for n{e) follows by setting C := '^'(0)- Next, if g E E^ is an arbitrarily fixed 
solution to (1 — Qo)g = y, then any other 77 e E^ with (1 — Qo)r] = y can be written uniquely in the form 
Va V ^ Q + '^B for some a S M. Writing w :— z'{Q) G Ei we have w ~ no(-, Q)ria — K^h by ( 12.201 ) 
and ( 12.111 ) with a e M determined by the constraint that w must belong to ker(L)^. This is obtained by 
observing that 

= Pfcw; = (t(no(-,0)f?) + a - k^{Koh))ll^{- , 0)B , 
that is, a = i<i{Koh) — ^(no(-, O)^'). For this a we put ^ := ija and get 

u{e)=eno{-,0)B + e^{nQ{-,0)^^Koh)+eu,{e) 
with |1w*(£)||ei = o(e2) as |£| ^ 0. 

□ 

2.3. Positive Solutions. We shall give conditions under which the nontrivial equilibrium solutions are 
positive. To this end we suppose that 

Eq is ordered by a closed convex cone Eq . (2.21) 

Then the interpolation spaces Ea- are given their natural order induced by the cone i?+ := E^ D E^ . For 
information on positive and strongly positive operators we refer to |j9]|25l. If {n, u) is a solution to i2.2i . 
(12.3b we say that u is a positive equilibrium provided that u{a) e Eq for a G J. 



Before turning to positive solutions we remark the following about the assumptions on Qq in Theo- 
rem |Z4l 



Remark 2.7. Assume that the parabolic evolution operator Iio{a, a) corresponding to Aq in ( 12.5b is pos- 
itive, that is. Ho (a, cr) G C+{Eq) for Q < a < a < am- If also G £+(Ei,£'tf) in ( 12.4b , then 
Qo G 1C+{E^) and thus the Krein-Rutman theorem entails that the spectral radius r(Qo) is (if nonzero) an 
eigenvalue of finite multiplicity with a positive eigenvector B G E^. Hence, in this case the assumption in 



Theorem \2.4\ that the normalized spectral radius r((5o) — 1 is an eigenvalue is no severe restriction. More 
restrictive is the assumption that this eigenvalue has geometric multiplicity 1. However, if Qq is strongly 
positive or irreducible, then t^Qq) = 1 has geometric multiplicity 1, see for instance \9, Sect. 12] and 1251 
App.3.2]. We also refer to the next section for concrete examples. 

Proposition 2.8. Suppose the assumptions of Theorem \2.4\ and ( I2.21I ). In addition, 
for each u E let A(u) generate a positive parabolic 

(2.22) 

evolution operator Ilu(a, a) , < a < a < a„j , on Eq ■ 
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If (n{e), u{e)) is the branch of solutions from Theorem \2.4\ then u{e) is positive provided that e £ (0, Eq) 
is such that 

-iu{e) = B + "fz{e) e E+ . (2.23) 
e 

In particular, if B belongs to the interior of , then u{e) is positive for e > sufficiently small. 

Proof. This follows from the fact that under the stated assumptions any solution (n, u) to (I2.2l i. ( 12.31 ) satis- 
fies 

u{a) ~ n„(a, O)u(O) , a ^ J , 

hence u{a) e E+ for a £ J if = w(0) £ E+ , and due to the fact that z{e) ^ in Ei ^ BUC{J, E^) 
as £ ^ 0. □ 

Remark 2.9. Recall that according to ^ II.Cor.4.4.2, II.Thm.6.4.2], A(u) generates a positive parabolic 

evolution operatorIlu{a, a) on Eq provided that A(u) £ C''{J, 7Y(£'i, Eq)) for some p > and — A(u)(a) 
is resolvent positive for each a ^ J. In this case, a solution u G Mi to (I2.2l i. (I2.3l l possesses additional 
regularity, see ^ II. Thm. 1.2.1, II.Thm.5.3.1]. 

Under some symmetry conditions on A and £ the equilibrium solutions provided by Theorem 12.41 are 
positive for each parameter value n{e), —Eq < e < Eq. More precisely, we have: 

Proposition 2.10. Suppose the assumptions of Theorem \2.4\ ( 12.21b , and ( I2.22l i. Let A{u) = A{—u) 
and £{u) — ~£{—u) for u £ E. Given m £ S sef QuW £{llu{-,0)w), w £ Ea, and suppose that 
Qu £ C,+ {Ea) for some a £ [0, Moreover, suppose that any positive eigenvalue of Qu has geometric 
multiplicity 1 and possesses a positive eigenvector Then 

C+ {{n{e),u{e)) ; ^u{e) £ E+] U {{n{e) , -u{e)) ; 7"(e) ^^^o""} 

consists of positive equilibria only. 

Proof Let e £ (— Eq, £q) \ {0}. Since u{e)) satisfies 

"(e) = n„(£)(-,0)7u(e) , 7u(e) = n(£)(5„(e)7w(£) , 

it follows that n{e)^^ > is an eigenvalue of Qu{e) with eigenvector 7m(£). By assumption there is a 
corresponding positive eigenvector Bu[e) and £ K \ {0} such that ju{e) — a^Bu{e)- If > then 
7u(£) £ Eq and thus u{e){a) £ E'^ for each a E J. Otherwise, if < 0, then —u{e) is a positive 
equilibrium solution with parameter value n{e) due to j{~u{e)) = —a^B^i^^-^ £ E'^ and owing to the 
symmetry conditions put on A and £. □ 

Proposition 12 . 8 1 guarantees that a branch of positive equilibria bifurcates from the branch of trivial equi- 
libria (n, u) = {n, 0), n £ R, at the critical value n = 1. Near the critical value n — 1 the set of n values 
corresponding to positive equilibria on the branch from Theorem 12. 41 consists of n values greater (i.e. su- 
percritical bifurcation) or less (i.e. subcritical bifurcation) than 1 depending on the sign of A(£) = n{e) — 1 
for £ > sufficiently small. If to > 2 in Theorem 12.41 this "direction of bifurcation", that is, the cases 
n{e) > 1 and n{e) < 1 for £ > small, depends on the sign of A'(0) = C (if nonzero), which in turn de- 
pends on (A^,(0))no(-, 0)B and Dl£^{0)[no{-, 0)B, Uo{-, 0)B] according to Proposition|22] Further, 
Proposition l2. lOl warrants under the symmetry conditions imposed that for any of the values n{e) / 1 there 
is a positive nontrivial equilibrium. Examples to which Propositions 12 . 8 1 and 12 . 1 Ol applv will be given in the 
next section. 



Under additional assumptions we can get more information about the positive equilibria and the direction 
of bifurcation. For simplicity we demonstrate this when £ is given by 

£{u) [ b{u,a)u{a)da , u £ E , (2.24) 
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where b e C""(I], L+(J)) with 1/p+l/p' = 1 and b{u, a) := b{u){a). Then (IZTSl l and ( |ZT91 ) clearly 
hold by putting 

€o(w) := / bo{a)u{a) da , i^,{u) :— / 0)7/(0) da (2.25) 

Jo Ja 

for 60(a) ■— b{Q,a) and 6*(ii,a) b{u,a) — bo{a). Let the assumptions of Proposition 12. 8 1 be satisfied 

and suppose that there exists e, G (0, Eq) such that (12.231 1 holds for e £ (0, e*). Let (I2.22l i hold and, given 

w e E, assume that 



Qu ■= / 6(u,a)n„(a,0)da , 
Jo 

belongs to JC+{E^). Note that Q^, for u — coincides with Qq defined in ( |2.10t . Set 

Ni := inf r(Q„) , := supr((5„) , 

where T := {m(£) ; e 6 [0, e,)}. Then 0<iVi<l<iVs<oo since r((3o) = 1. Moreover, 

nr{Qu)>l for {n,u) e A := {{n{e),u{e)) ; e £ [0,e,)} . (2.26) 
Indeed, given (n, it) G A \ {(1, 0)} we have u{a) — n„(a, O)u(O) for a £ J and 

^ u(0) = n£{u) = nQuu{0) , 
that is, 1/n is an eigenvalue of Q„ G C{E^), whence r{Qu) > 1/n. Suppose in addition that 

for each m G E , r{Qu) > is the only eigenvalue of Q„ G JC+{E^) with positive eigenvector . (2.27) 
This holds, e.g., if Qu is strongly positive. Then 

nr{Qu) = l, {n,u)eA. (2.28) 

Furthermore, letting 

[(Tj,o-s] := dM{n{e) ; e G [0,e,)} , 

it readily follows from (12.28b that 

0<cr, ==^<l<a, = ^<cx). (2.29) 

Therefore, under the assumptions of Proposition |2]8] ( 12.241 ). (I2.27l i. and if r((5„) < 1 for u G E, we have 
Ns < 1, hence 1 ^ Ng ^ ai and bifurcation must be supercritical in this case. Again, we refer to the next 
section for concrete examples. 

3. Applications to Population Dynamics 

We now apply the obtained results to problem ( |2.2| i, ( 12.31 ). Suppose ( |2.21| i and that the interior int(i?+) 
of i?+ is nonempty. Let A be of the form 

A(u, a) :— /i(u, a) + A{u, a) 

and 

£{u) := / b{u,a)u{a) da . 

Jo 

As observed in the previous section, £ satisfies (12.18b . ( 12. 19b with ( 12.25b provided that 

&gC"(E,L+(J)) , bo:^b{0)^0, (3.1) 

for some m > 1 and some ball E in Ei centered at with radius Rq > 0. Moreover, regarding Proposi- 
tion |2j6] we note that £'„£*(0) = 0. Let a G [0, <;) and let $ be the ball in Ea with center and radius 
R>0. Let 

AeC"\^,C{Ei,Eo)) (3.2) 
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be such that 

— A{w) generates an analytic semigroup on Eq and is resolvent positive for each w e $ . (3.3) 

Making i?o > smaller if necessary it follows from the compact embedding Ea and ( 12. Il l analo- 

gously to ill] Thm.6.2,Thm.6.4] that the Nemitskii operator of A (again labeled A), given by 

A{u){a) := A{u{a)) , a^J, u G S , 

belongs to C™(E, Loo( J, C{Ei,Eo))). Since Ei ^ BUC''-^{J, Es) for 5 e [0, <;) owing to (EB and the 
interpolation inequality |l2] I.Thm.2.1 1.1], we deduce from (13.31 ) and Remark ing] that 

[a^ A{u{a))]^C'-''{J,n{Ei-E^)) 

generates a positive parabolic evolution operator Ua{u){0', o") on Eq for each u G S. Set Aq :— A{0) and 
suppose there exist llJq > and cf) G [0, 7r/2) such that uiq > type(— Aq) and 

ujo + Ao£BJP{Eo;cj)). (3.4) 

Moreover, suppose that 

g-aAo g strongly positive for a > . (3.5) 

If /i is a function such that 

[^^^(u,.)] gC"(E,L+(J)), (3.6) 
we set fj,o{a) :— fi{0, a) for a G J and further suppose that 

HO e BUC{J) , inf ^o(a) > (3.7) 

aG J 

and 

lim /io(a) exists if a„i = oo . (3.8) 

a — ''Oo 

Put A(m, a) := a) + A{u{a)) for a G J, u G S and note that 

Ao(a) := A(0, a) = iio{a) + , a £ J . 
Clearly, A{u, •) generates a positive parabolic evolution operator n„(a, a) on i?o for each u G E given by 

n„(a, d) := e- {7^(„) {a,a) , < a < a < a,„ . 

From (|14] |. dlTT i. and fT III.Cor.4.8.6] it follows that we may apply Remark|23]b) to conclude that ( 12.151 1 
holds true. Proposition l2.2l now guarantees that A satisfies (12.5b provided Eq is an UMD space. 
Finally, let Qq G IC+{E^) be given by 







and note that ( 13. 5t and feo ^ imply that Qq G IC{E^) is strongly positive, hence irreducible (see ||9] 
Sect. 12]). In particular, since the interior of E^^ is assumed to be nonempty, it follows from |9, Thm.12.3] 
that r{Qo) > is a simple eigenvalue of Qo with a corresponding eigenvector B G int(i?+) (and this is the 
only eigenvalue with a positive eigenvector). Let then 6o be normalized such that r{Qo) = 1. 



Combining Proposition 12. 21 Remark l23l Proposition l2.8l and Theorem l2.4l we obtain: 

Theorem 3.1. Let Eq be a UMD space satisfying (I2.21l i and let int{E^) ^ 0. Suppose (I3.1l l-( l3.8l ). Then 
the problem 

daU + A{u{a))u + a)u = , a G J , 



u{0) — n / b{u, a)u{a) da , 
Jo 

has a branch of nontrivial solutions (n(£),u(e)) G K+ x Ei, < |e| < So, of the form 

n{e) = 1 + A(e) , u{e) = e(no(-, 0)B + z{e)) 
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such that A : (— eoi £o) ^ M, z : (— £oj £o) ^ lEi are m-times continuously dijferentiable and A(0) — 0, 
z{Q) = 0. If e > w sufficiently small, then u{e){a) G for a E J. 

If, in addition, the symmetry conditions 

A{~u) = A{u), ^ii-u, ■) = fi{u, ■) , b{~u, ■) ^ b{u, ■) (3.9) 

hold for u € T, and if 

ram 

b{u, a) e- ^0 M(u,r)dr ^^^^^ 

for u G S is such that Qu G C+{E^) and 

any positive eigenvalue of has geometric multiplicity 1 

and possesses a corresponding positive eigenvector , 
then it follows from Proposition |2.10| 

Corollary 3.2. Suppose the assumptions ofTheorem \3.1\ tosether with ( 13.91 ), ( 13.1 Oi l. Then, for each param- 
eter value n(e), e G (— eoi ^o) \ {0} provided by Theorem \3. 1 1 there exists a positive nontrivial equilibrium 
solution of the form u{e) or —u{e). 

3.1. Example. Let Q, C M^, > 1, be a bounded and smooth domain lying locally on one side of dVl. 
Let dil = To U Fi, where Fq, Fi are both open and closed in dft and Fq n Fi = 0. Consider 

A{u, x)w :— — • [a{u, x)Vxw) + ai{u, x) ■ WxW + oo(m, x)w , 

where 

[u^a{u,-)] G C"($,Cl+"(n)) , 

for some to > 1, G (0, 1) small, and some open ball $ in C^^°'{ft) around 0. Moreover, assume that 

a{u, x) > , x efl , u e ^ . (3.12) 

Let 

ryoeC^Ti) (3.13) 

and let denote the outward unit normal to Fi. Let 

w , on Fo , 

+ vo{x)w , on Fi . 

Fix p > iV + 2 and let i?o := Lp{fl) he ordered by its positive cone of functions that are nonnegative almost 
everywhere. Note that Eq is a UMD-space. Set Ei g(ri), where 

r w^^[n) , o<2e< i/p, 

-^l {w(^ W^^iVL) ■ M|r„ = 0} , 1/p < 2^ < 1 + 1/p , 
[ [w G W^^{n) ; = 0} , 2f > 1 + 1/p , 

are subspaces of the usual Sobolev-Slobodeckii spaces Wp^{Vl). Note that if E'^ :— {Eq, Ei)^^p, then 

= K'b ^E^= ^C^+^i^), l + n/p + a <2a<2q , 
where dots indicate equivalent norms l|27l . and int(£'+) ^ 0. Consider 

A{u)w -.^ A{u,-)w , w & El , uG$. 
Then i3.2\ and (I3.3l l follow, e.g., from ([T]. Moreover, suppose that 

flo(0,2:)>0, ai(0,a;)=0, xgH, 

/ ^ in (3.14) 

iyo{x) > , X G 1 1 , 



B{x)w 
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and put Aq A{0). According to 1 1 1, —Aq is resolvent positive and generates a contraction semigroup 
on each Lq{il.), 1 < q < oo, is self-adjoint in £2(^1), and there exists a largest eigenvalue Ao < of 
—Aq € C{Ei, Eq) with a positive eigenfunction B E . Moreover, jO] Cor.13.6] ensures (13.51 ). From 
||2l III.Ex.4.7.3.d)] we deduce now (|3.4| i for each o^o > 0. Given G (0, 00] and some open ball E in 

El = Lp(J, W^gl^)) n W^i J, Lp{n)) centered suppose that 

satisfies , <E3 with luq = , and ■ (3.15) 
Thus, if bo '■= b{0, •) for b := [u ^ b{u, •)] e C""(S, L^, (J)) is nontrivial and normalized such that 

/ "&o(a)e-/o>o(r)drgAoaj^^ -^^ (3^^^^ 

Jq 

then e'^^'^B = e'^^°B entails QqB = B, where 

Thus r((5o) = 1 by ( |3.16t since r{Qo) is the only eigenvalue with positive eigenfunction. Therefore, 
Theorem |3.1| entails : 

Proposition 3.3. Let p > N + 2 and suppose (13.11) , ( I3.11b -( l3.16b . Then the problem 

daU + A{u{a)^ x)u + ^{u, a)u = 0, a £ J , xefi, 

u(0) = n b{u, a)u{a) da , x £ , 
Jo 

B{x)u = 0, a>0, X edn, 
has a branch of nontrivial solutions 

(n(e),w(e)) e M+ x (Lp(J, n J, Lp(f]))) , < |£| < Eq , 

bifurcating from (n, u) = (1, 0), iMc/i that u{e){a) G L^(fl) for a £ J and e > small. 

Remark 3.4. The proposition above also holds if Eq :— Lq{^) and Ei := for q > N + 2 

different from p £ (1, oo). The only difference is that the interpolation space E^ equals a subspace of the 
Besov space Bq''p{fl). 

3.2. Example. We may also consider a functional dependence of A on u. Indeed, let again G (0, oo] 
and let fi, Ei, and Eq be as in the previous example with p G (l,oo) arbitrary. Given w G Ei 
Li(M+, Lp(f7)) let U u{a)da and consider aIu)w A(U, ■)w for w £ Ei ^ with A, B 

as in the previous example satisfying ( |3.12| )-( [3T4l i but $ in (13. lit is now an open ball in Lp{il) centered at 
0. Suppose dTTSl l with infa>o a) > type(-A(w)) for w G S, dTBI l, (lITel l and 6 G C™(I], L+ (J)). 
Moreover, assume that b{u) ^ 0. Then, analogously to the previous example, 

Qu := / " 6(M,a)e--''o°''("''-)'i'-e~'^'^("Ma G /C+(£;,) 

is strongly positive for each w G S, whence (I2.27l i by |j9] Thm.12.3, Cor.13.6]. We obtain from Theorem l3 . 1 1 
a branch of nontrivial solutions 

{n{e),u{£)) GR+ X {LpiJ^WlsiQ)) nW^iJ, Lp{n))) , < |e| < , 

to the problem 

daU + A{U, x)u + a)u = 0, a G J , x £ il ^ 

u(0) — n b(u, a)u{a) da , x G 17 , 
Jo 

B{x)u = 0, a>0, xGarj, 
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bifurcating from (rt, u) = (1, 0), such that u{e) is positive for e > sufficiently small. If Xq{u) denotes the 
largest eigenvalue of —A{u) E C{Ei , Eq) for u G S and if 

'6(u,a)e^^°^("''')'^'^e"^"(")da < 1 , ugE, (3.17) 

Jo 
then 

r{Qu) < 1 , M e S . 

Indeed, if i?„ is a positive eigenfunction corresponding to Ao(u), then e^°'^'^^^Bu — e^'^'^^'^Bu entails 

g„B„ = biu, a) e- ^0 e^"^")" da , 



whence r{Qu) equals the left hand side of (I3.17l i. Recalling (|2.29l l we deduce that bifurcation must be 
supercritical provided ( I3.17l i holds; that is, for e > small we have n{e) > 1 and u{e) is nonnegative. Note 
that Ao(u) < if ao{u, •) > and ao{u, •) — div(ai(M, •)) > in 17, z^o > and ai(u, •) • > on Fi (see 
[l , Rem. 11.3]) in which case the term e'^o(f)a jjj (|3,17| | can be neglected. Moreover, type(— < 
in this case. If the functions a, ai, ap as well as fi and b are symmetric with respect to u, that is, if 
a{u, •) = a{—u, •) etc., then Proposition 12. 1 01 entails that there is a positive equilibrium solution for any 
value of n{e), — £o < e < £o- 

3.3. Example. Let G (0, oo] and let f2 C M^, > 1, be a bounded and smooth domain. Given 
w G El let U :~ "* u{a)da and consider 

A{u)w -V, • («([/, x)\7,w) , w e Ei -.^ W^^in) {v G W^in) ; d,v = 0} 

so that -A{u) G H{Ei,Eq) for := Lp{n) andp G (l,oo) provided that a G C"' {Lp{n) , (0.)) is 
such that a{U, x) > for x E ft. If {n, u) is any positive solution to the problem 

daU + A{u)u + a)u = , a > , 



u{0) = n / b{u, a)u{a) da , 
Jo 

then the relation u{0) = nQuu{Q) must hold, where 



Therefore, owing to the fact that 

e-'^'^("Vdx / (f>dx, (j)eLp{n), (3.18) 



n 



it follows by integrating the previous relation that necessarily 

l=n f " a) e--/'o°^("^'')'^'' da (3.19) 



which is the same constraint as in the non-diffusive case (see f6\). Let b G C™(Ei, L^,{J)) and suppose 
(|3.15t . (13.161 1 with Aq = (in particular q{0) ~ 1). Further assume that 

b{u, a) < 6(0, a) = &o(a) , /^(w, a) > /i(0, a) = /io(a) (3.20) 

for a E J and u G E^^, which is a common modeling assumption stating that effects of population densities 
do neither increase fertility nor decrease mortality. Then q{u{e)) < q{0) — 1 for the positive solution 
{n{e),u{£)), £ > small, provided by Theorem 13. II Thus ( |3.19t entails n{e) > 1 for e > small, 
that is, bifurcation must be supercritical, and there is no equilibrium solution other than the trivial u = 
corresponding to a parameter value n < 1. 
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We shall point out that the present example simply reflects the non-diffusive case in the sense that our 
results here actually follow from the case A = (see 161). For this it is enough to observe that Aq = is an 
eigenvalue of —A{u) with corresponding constant eigenfunctions. 

Moreover, taking i? = 1 we have Ilo{a,0)B — e^-^o f^o{s)ds ^j^^ since the projection onto M-^ in 
Proposition 12.61 is given by 

1 - Pm = ■^(^) ' 

the direction of bifurcation, given by ( in Proposition 12.61 can in principle be computed explicitly using 
(13.181 1 also if one does not assume ( |3.20| i. 

Remark 3.5. In all our examples we omitted a dependence of ii and b on the spatial variable for simplicity. 
However, it is clear that such a dependence can be included as well. 
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